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Abstract. We analyze the joint extremal behavior of n random products of the form Yl'jLi X"' j , 1 < i < n, for 
non-negative, independent regularly varying random variables Xi ,..., X m and general coefficients (Hj £ R. Products 
of this form appear for example if one observes a linear time series with gamma type innovations at n points in time. 
We combine arguments of linear optimization and a generalized concept of regular variation on cones to show that 
the asymptotic behavior of joint exceedance probabilities of these products is determined by the solution of a linear 
program related to the matrix A = (aij). 
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1. Introduction 

The tail behavior of products of powers of heavy-tailed positive random variables is of crucial importance 
in many applications, particularly in finance, but e.g. in network modeling too. In stochastic volatility time 
series, the log-volatilities are usually modelled as linear time series 

OO 

log crt = ^2 a iVt-i , te z, 

2—0 

If the innovations rji,i £ Z, have an exponential or gamma type tail, then the volatility a t at time t is a 
product of powers of the regularly varying random variables Xi := e Vi , i £ Z, with exponents depending on t. 
To assess the risk of a volatile market at different time points t\,...,t n , one thus has to analyze probabilities 
of the type -Pdlyli AA > x > 1 < i < n) for suitable exponents a^. Using a new Breiman type result, it 
was shown in Janfien and Drees (2016) that these probabilities also determine the risk of jointly large losses 
over different periods. 

Similarly, in a credit risk model for n risks with k independent factors Z \,..., Zk, the i-th risk is often 
modeled as a multiple of exp(Xy = i ®ijZj + T)) with Y),l < i < n, denoting the idiosyncratic part (cf. 
Embrechts et al. (2014)). If the Z 3 , 1 < j < k, and Yi, 1 < i < n, have an exponential or gamma type tail, 
the analysis of the joint tail risk again leads to probabilities of the above type. 

In network modeling, both transmission durations L and rates R arising from one source may be modeled 
by regularly varying random variables with different indices az, and an (see, e.g., Maulik et al. (2002)). The 
total volume of traffic from one source can then be expressed as X^X^ * 1 for random variables 
which are regularly varying with index — 1. If one wants to determine the probability that different sources 
contribute large volumes in the same period, then again probabilities of the above type arise. Moreover, one 
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may introduce dependencies between Xl and Xu for the same or for different sources by modeling them as 
products of (partially) identical factors with different exponents. 

As the example of log-volatilities demonstrates, the analysis of the joint tail behavior of power products is 
equivalent to the corresponding analysis for linear combinations of random variables with exponential type 
tails such that exponentials of these random variables are regularly varying. Hence the results given below 
allow a tail analysis in such settings too. 

Motivated by these examples, we analyze the asymptotic behavior of the probability of joint exceedances 
of power products, i.e. of 

> CiX, 1 < i < , Ci > 0,1 < i < n, (1-1) 

where Xj are independent, non-negative regularly varying random variables and a.y are coefficients which 
may be negative. We restrict the analysis to the case of a finite number m of factors, but extensions to 
an infinite number of factors, using a Breiman-type argument, are possible. In Janfien and Drees (2016), 
such probabilities were investigated under the restrictive assumption that no coefficient is negative. It was 
shown there, that the probabilities behaved asymptotically like a multiple of JI, P(X, > x Ki ), where 
k = (k\, ..., n m ) is the solution to a linear program determined by the matrix A = (o-ij). While the 
restriction to positive coefficients a, t j seems acceptable e.g. for multi-factor models, it is quite severe for log- 
volatility time series. Moreover, essentially only the case n = 2 was considered in Janfien and Drees (2016) 
and the techniques employed do not easily generalize to higher dimensions, which limits the applicability of 
the established results further. 

Using a recently introduced abstract concept of regular variation on cones based on the notion of M- 
convergence (see Hult and Lindskog (2006), Lindskog et al. (2014)), we can avoid all these drawbacks. 
To this end, we first introduce a non-standard form of regular variation on the cone (0,oo) m for the 
random vector (X \,..., X m ), from which one may conclude the asymptotics of probabilities of the type 
P((X i/x Kl ,... ,X m /x Km ) £ B) for suitable coefficients Ki,...,K m and sets B C (0, oo) m that are bounded 
away from the boundary of the cone. Unfortunately, in general the sets M = {x £ x °j') > a, 1 < 

i < n} pertaining to the probabilities (1.1) are not of this type. Hence, quite involved arguments are needed 
to prove that the parts of M close to the boundary of the cone are asymptotically negligible. To this end, 
auxiliary results are proved in Section 4.3 which are of interest on their own. In particular, Proposition 4.4 
can be seen as a multivariate version of the direct half of Karamata’s Theorem, cf. Remark 4.5. 

The outline of this paper is as follows: In Section 2 an abstract notion of regular variation on cones is 
briefly recalled, with a view towards the later application of this concept. Our main results are stated in 
Section 3, with Theorem 3.3 being the central conclusion. Proofs of the results are given in Section 4.2 while 
some auxiliary results needed in the proofs are gathered in Section 4.3. 

Notations and conventions 

We write bold letters for vectors, i.e. x is short for (aq,..., x n ) £ ffi" if it is clear that x is of dimension n. 
The i -th component of x is denoted by aq. We write 0 and 1 for a (column) vector of suitable dimension 
which consists of only zeros or only ones. Inequalities for vectors are meant to hold componentwise. We 
denote the complement of a set A by A c and its boundary by dA. For x £ I™ and A C K” we set d(x, A) = 
inf a6j 4 ||x — a||, where || • || denotes the Euclidean norm. Similarly, we set d(A,B) = infag^bes || a — b||. For 
A C K" and r > 0 set A r := {x £ ffi" : d(x, A) < r}. Denote the Borel sigma algebra on R" by B” and for 
a set A G B n write B n fl A = {B £ B n : B C A}. We write A(-) for the Lebesgue measure on B n . 

2. General regular variation on cones 

In the following, we will make frequent use of an extension of the concept of multivariate regular variation 
which was introduced in Hult and Lindskog (2006) and Lindskog et al. (2014). For some m G N, let ® : 
(0, oo) x [0, oo) m [0, oo) m , (A, x) >->■ A®xbea “multiplicatish” mapping with the following two properties: 
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(Al) the mapping ® is continuous, 

(A2) 1 <g> x = x and for Ai, A 2 > 0 we have Ai <g» (A 2 ® x) = (Ai ■ A 2 ) <g> x for all x £ [0, oo) m . 

Consider a closed subcone C of [0, oo) m w.r.t. this mapping, that is, A <8> C := {A <8> x : x £ C} C C for all 
A > 0. We assume that the following condition holds: 

(A3) d(x, C) < d( A ® x, C) if A > 1 and x £ O. 

The complement O := [0, oo) m \ C is an open cone, which is assumed not to be empty. 

The notion of regular variation on O w.r.t. <8>, which is introduced below, rests on the definition of 
convergence in the space Mo of Borel measures on (O, B m (~l O) whose restrictions to [0, oo) m \ C r are finite 
for each r > 0. Denote by C + (0) the class of non-negative, bounded and continuous functions / on O 
vanishing on C r for some r > 0. We endow Mo with the topology that is generated by open sets of the form 


v £ 


fidv - / fid/i 


< e, 1 < i < k 


with n £ Mo,/,; £ C + (0),i = 1 and e > 0. A Portmanteau Theorem (cf. Lindskog et al. (2014), 

Theorem 2.1) shows that convergence of measures v n to a measure v in this topology is equivalent to 
the convergence v n (A) —> v(A) for all Borel sets A in O which are bounded away from C and for which 
v{dA) = 0. 

Definition 2.1 (see Lindskog et al. (2014), Definitions 3.1 and 3.2). A measure v £ Mo is called regularly 
varying on O with respect to the mapping ® if there exists an increasing, regularly varying function c : 
[0, 00 ) —>• (0, 00 ) and a nonzero measure ft £ Mo such that 

c{x)v{x <g) •) —> p(-) in Mo as x —> 00 . 

Lemma and Definition 2.2 (see Lindskog et al. (2014), Theorem 3.1). Definition 2.1 implies that there 
exists an a > 0 such that 

fj(X <g> A) = X~ a p,(A) (2.1) 

for all A > 0 and Borel sets A C O. We call —a the index of regular variation of the measure v in Definition 
2.1. The value of a in (2.1) is equal to the index of regular variation of the normalizing function c in Definition 
2 . 1 . 


Proof. Equation (2.1) is stated in Lindskog et al. (2014), Theorem 3.1. By this and (A2), we have, for all 
A > 0 and A in O which are bounded away from C and for which v(dA) = 0, that 

c(Ax) v(x®A) v(x® A) p(A) 

Inn , , = lim . -— = Inn —-—-— = ——-— = A . 

x-s-oo c(x) x->oo i>([Ax) <g) A) x->oo v(x igi (A ® A)) /i(A <g> A) 

Therefore, c is (univariate) regularly varying with index a. □ 

Definition 2.1 unifies several different concepts of regular variation of a random vector X = (Xi ,..., X m ) 
with values in [0, oo) m and distribution v. 

Example 2.3 (Multivariate regular variation). If ® denotes the usual scalar multiplication A®x = Ax and 
C := { 0 }, then O = [0, oo) m \ { 0 } and Definition 2.1 reads as 

as x —> 00 for all p-continuity sets A C B m n [0,oo) m bounded away from 0. This is the classical regular 
variation ofX (see e.g. Resnick (2007), Section 6.1.f). 
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Example 2.4 (Ledford-Tawn-model). If C := [0,oo) m \ (0, oo) m , then regular variation on O = (0,oo) m 
w.r.t. the usual scalar multiplication has been considered by Ledford and Tawn (1997) in the bivariate case 
m = 2 (after suitable marginal standardization). It is equivalent to the convergence 

c(x)P ^ ..., £ Aj -> p(A) 

as x —t oo for all p-continuity sets A C B m PI [0, oo) m bounded away from both axes in the case m = 2 resp. 
from {x : Xj = 0 for some 1 < i < m} in the general case. 

Note that a random vector X may be regularly varying in the classical sense of Example 2.3 and in the 
present sense with different normalizing functions c resp. c. If c(x) = o(c(x)) as x ^ oo, then X is said to 
exhibit hidden regular variation (cf. Resnick (2007), Section 9.).1). 

Here, we consider a different mapping <g> and different cones as well. Let, for k £ [0,oo) m , 

: (0, oo) x [0, oo) m -?• [0, oo) m , 

(A, (xi , • ■ - , Xm) ) 1 t A (xi, ... , Xm ) * = (A 1 X\ ,..., A m X m ). 

We want to analyze the asymptotic behavior of m-dimensional non-negative random vectors that have ex¬ 
treme values inn 6 {1,..., m} of their components. For ease of notation, assume that the first n components 
of n are positive and the last m — n components are equal to zero, so that x Ki —^ oo as x —> oo only for 1 < 
i <n. Define the cones C n = ([0, oo) n \ (0, oo) n ) x [0, oo) m ~ n and O ra = [0, oo) m \C n = (0, oo) n x [0, oo) m ~ n 
w.r.t. the mapping ® K . Since 

d(x, C n ) = min{xi,..., x n } < min{A K1 Xi,..., A Kn x„} = d (A <g> K x, C„) 

for all A > 1 and x £ 0„, the assumptions (Al)-(A3) are satisfied. Note that in the case n = m and n = 1, 
the regular variation on O n w.r.t. <S) K is equivalent to the concept of regular variation considered in Example 
2.4. 

Lemma 2.5. Let k. £ [0, oo) m with Ki > 0,1 < i < n, and Ki = 0, n < i < m, for some n < m. Furthermore, 
let X\,... ,X m be independent, non-negative random variables such that X\,... ,X n are regularly varying 
with index —1. Then, p( x Ai<i<m regularly varying on O n w.r.t. <S> K with index —a = — Ki- 

Proof. From the independence of Xi, ..., X m and the regular variation of X\, ..., X n it follows that 


lim 


= lim 

x—foo 


p(Xi) i<j<m I X <g) re (I X { a ii oo) I x I X [bj, oo) 

\ \ / \n<j<m y 


rnu p(Xi>x^) 

nr=i nxi > aix Ki ) riyLn+i 


nr=i p(x i >x«) 


= H*7 1 n P(Xj>bj)=:0 

2=1 j=n -\-1 


X (di,oo) I x I X [bj, oo) 
^l< 2 <n J XTiKj^m / 


( 2 . 2 ) 


for all a, > 0,1 < i < n, and bj > 0 ,n < j < to. Since these limits are finite, we have shown that the 
family of measures c(x)P( a -j)i<)<>»» (x <8> K ■), x > 0, is relatively compact in Mq„ (cf. Theorems 2.4 and 2.5 
in Lindskog et al. (2014)), where c(x) = (Jl" =1 -PpQ > x Ki )) '. Furthermore, all accumulation points of 
this family agree on a generating 7r-system. Thus, p( x Ai<j<m j s regularly varying on O n w.r.t. The 
index of regular variation follows from Lemma and Definition 2.2, since c is regularly varying with index 

V" k — V m k n 
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3. Joint extremal behavior of random power products 

In the following, let X\,..., X m ,m £ N, be independent, non-negative random variables, not necessarily 
with the same distribution. We will give asymptotics for the joint exceedance probabilities (1.1) of n < m 
“power products” 

m 

n^r, 1 <i<n, (3.1) 

3 =1 

over the same threshold x as x —> oo for rather general values of ajj £ R, 1 < i < n, 1 < j < m. A product 
may take the value +oo if Xj = 0 and atj < 0 for some i. j, but throughout we use the convention that 
+00 -0 = 0 and 0° = 1. 

In order to derive our results we make some assumptions about the tail behavior of the Xj, 1 < j < m. We 
assume that all or at least the “relevant” (in a sense specified below) Xj, 1 < j < m, are regularly varying 
with index —1. We will see that the joint extremal behavior of the products in (3.1) is closely related to the 
solution of the linear optimization problem 

m 

find x > 0 such that Ax > 1, E Xi —>■ min! (3-2) 

i= 1 

with A = (aij)i<i<n,i<j<m £ M nxm . Before we give proofs for the asymptotic behavior of the joint ex¬ 
ceedances, we want to motivate the connection of our question to the linear optimization problem in (3.2). 
To this end, assume for simplicity that > 0 for all 1 < i < n, 1 < j < m. Let y > 0 be a feasible solution 
to (3.2), i.e. Ay > 1. Note that for all x > 1 

m 

Xj > x Vj , 1 < j < m =>- Xj 1,3 > x^i =1 aijVj > x, 1 < i < n, 

3 =1 


by (3.2) and thus 


P 


Xj' 3 > x, 1 < * < 


> n p ( x > > * w )- 

3 = 1 


(3.3) 


If all Xj, 1 < j < m, are independent and regularly varying with index —1, the right hand side is a regularly 
varying function in x, with index a = — Vj- Now, the smaller the value of |a|, the slower is the decay 

of the function on the right hand side as x —>• oo. So, heuristically, if the value of |a| and thus the value of 
Y'J—i Vj is minimized, this is the most likely combination of extremal events for the single Xj which leads 
to joint extremal behavior of the power products (3.1). We will see in Theorem 3.3 that the right hand side 
of (3.3) is not only a lower bound for the joint exceedance probabilities but also, under some additional 
assumptions about real valued A, tail equivalent to it. For a general (not necessarily non-negative) matrix A, 
the next theorem gives upper and lower bounds for the order of decay of the joint exceedance probabilities. 

Theorem 3.1. Let X\,... ,X m be independent non-negative random variables. Let k = (sq ,...,K m ) T be 
an optimal solution to (3.2). 

(a) Assume that all Xj, 1 < j <m, are regularly varying with index -1. Then for all e > 0, 




= o 



^11 x T j > x ' 1 ^ ^ n j j > 


x oo. 


(3.4) 
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(b) Assume that E(Xj s ) < oo, 1 < j < m, for all S £ (0,1), and additionally that there exists c > 0 such 
that P(Xj > c) = 1 for all 1 < j <in with Kj = 0. Then for all e > 0, 

Xj' 3 > x, 1 < i < = o(x~^ i=lKi+e ), x —> oo. (3.5) 

Remark 3.2. In contrast to the other results of this and the following section, the assumptions of Theorem 
3.1 (b) do not include regular variation of at least some of the Xj, 1 < j < m. Note, however, that regular 
variation with index —1 of X ±,..., X m implies that E(Xj~ s ) < oo, 1 < j <m, for all S £ (0,1). 

The proof is given in Section 4.1. Under some additional assumptions about the structure of A, the 
following Theorem 3.3 gives precise asymptotics for the joint exceedance probabilities of the random power 
products. 

Theorem 3.3. Let A = ( aij ) £ R nxm ,n < m, be such that the optimal solution n to the linear optimization 
problem (3.2) is unique and non-degenerate (i.e. it has n positive components) and denote by A K £ R raxn 
the matrix which is derived from A by deleting all columns 1 < j < m for which Kj = 0. Then this matrix 
is invertible. 

Let X\,... ,X m be independent non-negative random variables and assume that there exists e > 0 such 
that for 



1 < j < 77i with Kj > 0 : Xj is regularly varying with index — 1 


1 < j < tu with Kj = 0 : 


E ^xj 1 Are A -b +e ^ < oo and 


Then 


lim 

x —^OO 


£l(xj lTA “ lA)i_6 ) < oo 
p (n™ix; y > c,x, i<i<rij 


n p(Xj>x«>) 

{l<j<Tn:Kj > 0 } 


- l “ 4 ' rl o^ 1 £ e ( A ' 


j:K j=0 


=: M ( X (Ci,oo) 

\l<i<n ) 


(3.6) 


(3.7) 


(3.8) 


for all Ci > 0,1 < * < n. In particular, the distribution of (nr=i*T) i <i<n is regularly varying on (0, oo) n 
w.r.t. scalar multiplication (in the sense of Example 2-4). The normalizing function can be chosen as c(x) = 
(n { i< 7 -< m . K .> 0 }^(*,' > a: Kj )) _1 and the corresponding limit measure is g, as above. The index of regular 
variation is equal to — X7=i = ~^ 

Remark 3.A. Given the expression in (3.8), it is obviously necessary to assume that E (xj 1 Are A ^ < oo for 
all 1 < j < m with Kj = 0 in order to ensure that the limit in (3.7) is finite. The actual assumption about the 
moments of those Xj with k 3 = 0 which is stated in Theorem 3.3 is similar to the assumption in Breiman’s 

lemma (cf. Breiman (1965)) in that we need a little more than the finiteness of the moments E ^xj 1 Are A ^ ^ 

in order to apply a dominated convergence theorem. Furthermore, note that for (l 7 A~ 1 A )j > 0 only the 
first moment assumption is necessary, the second follows for e > 0 chosen small enough. On the other hand, 
if (l T A f : 1 A)j < 0, then only the second assumption is necessary. Only in the case (l 7 A k . 1 A)j = 0 both 
assumptions are necessary. 

Remark 3.5. (a) Under the given assumptions the vector l 7 A^ 1 which appears in the statement of The¬ 

orem 3.3 is the transposed of the unique optimal solution to the so-called dual problem to (3.2), which 
is given by 


m 

find x > 0 such that A 7 x < 1, Xi —> max!, (3.9) 

i= 1 
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see the proof of Theorem 3.3 and (4.18) for details. This implies that ( 1 7 A^ 1 A) T = A 7 (l 7 A^ 1 ) 7, < 1 
and by the assumed uniqueness and non-degeneracy of the solution even (l 7 A f : 1 A) J - < 1 for those j 
with Kj = 0, cf. the remark after (4.27) in the proof of Theorem 3.3. Therefore, the assumptions in 
(3.6) are always satisfied if all Xj , 1 < j < m, are regularly varying with index —1 and bounded away 
from 0 (or their distributions concentrate sufficiently little mass around 0). 

(b) For a general linear program of the form 


m 

find x > 0 such that Ax > 1, ^f^ZjXj —> min! (3.10) 

l=i 


with optimal solution k, the value of Zj — (l 7 A~ 1 A)j is sometimes called the reduced cost of variable 
1 < j < m. If Kj = 0 in the optimal solution, then this solution is not affected by a change of Zj in the 
objective function in (3.10) as long as Zj > (l 7 A f T 1 A)j, cf. Section 2.5.1 in Sierksma (1996). In the 
context of Theorem 3.3, the values of (l 7 A^A)^ for j with Kj = 0 can be interpreted in a similar 
way, since the left or right tail behavior of these Xj does not influence the extremal behavior of the 
random products (except for a possible change in the multiplicative constant of the limit) as long as 
there exists e > 0 such that 

P(Xj >x) = 0{ ara^A^Ab-c), if (1 t A“ 1 A)j > 0 

PiX- 1 >x) = 0{x^ TA - lA ^~ e ), if {l T A~ 1 A) j < 0 ' 


The proof of Theorem 3.3 is given in detail in Section 4.2, but we want to lay out briefly the main idea 
here. To this end, assume w.l.o.g. that the first n components of k are positive and the last m—n components 
are equal to zero. The main idea is to use the regular variation of the measure v := <S>"=i P Xi (where P Xi 
stands for the law of Xi ) on O n w.r.t. <8> K , cf. Lemma 2.5. Furthermore, we show that under our assumptions 
the equality Ak = 1 holds. Then, 


P I X '“’ J > dx, 1 < i < n 


Vj =1 


m / -r r \ a. 


= P (n(^) !<*<») =v(x®kM) 


for 


M = M( A, c) := < x G [0, oo) m : x' A} > d, 1 < i < n > . 


j=i 


(3.11) 


The next step is to apply the Portmanteau Theorem (cf. Theorem 2.1 in Lindskog et al. (2014)) to show 
convergence of the right hand side in (3.11) under suitable normalization as x —> oo. Note, however, that 
the set M is not bounded away from C n (cf. Section 2), so we cannot directly apply this argument. As an 
intermediate step, we therefore have to show that we can replace M by M (~l (6, oo) n x [0, oo) m-n , <5 > 0, in 
(3.11) and that under the necessary normalization the difference is negligible as <5 \ 0. 

The following example illustrates the statements of Theorem 3.3 and in particular the role of negative 
coefficients a^. It also demonstrates applications of Theorem 3.3 for the extreme value analysis of time 
series. 

Example 3.6. Let be a log-linear time series of the form 


ln(y t ) =ln(X t )-0.51n(X t+1 ), t g Z, 


where X t ,t £ Z, are i.i.d. and regularly varying with index —1. Using Theorem 3.3, we can derive the 
asymptotics for the probability that three consecutive extreme observations of similar magnitude occur, i.e. 
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for P(Y l > c\X,Y 2 > c 2 x,Y 3 > c 3 x), Ci > 0,i = 1,2,3. Rewrite this probability as 

P(X iX-°' 5 > c lX ,X 2 X -° 5 > c 2 x,X 3 X 4 0 - 5 > c 3 x) = P j X“ ij > Ci x, l<i<3 

V =1 

with 

/ 1 -0.5 0 0 \ 

A = )i<i<3,i<j<4 = 0 1 —0.5 0 

\ 0 0 1 -0.5 / 

The optimal solution to (3.2) is then given by k = (7/4,3/2,1,0) and this solution is unique and non¬ 
degenerate. Furthermore, 



{ I 

' 1 

-0.5 

0 \ 1 \ 

(1 T A k 1 A) 4 = 

w 

0 

1 

-0.5 A = ((1,3/2, 7/4)A) 4 = —7/8 


1 

V 0 

0 

1 / A 


Let us first additionally assume that E(X 4 7 ^ 8 e ) < oo for some e > 0. Then all assumptions of Theorem 
3.3 are satisfied and the random vector {Yi,Y 2 ,Y 3 ) is regularly varying on (0,oo) 3 with respect to scalar 
multiplication. The index of regular variation is equal to — Sj=i K j = —17/4 and the limit measure p is 
given by 

M x Li(ci,oo)) = | r cr 1 c 2 - 3/2 c 3 - 7/4 i;(x 4 - 7/8 ). 

Note that the negative exponents do influence the solution of the optimization problem and hence the index 
of regular variation. If, for instance, in the matrix A —0.5 is replaced with —0.25 everywhere, the optimal 
solution is given by (21/16,5/4,1,0) and the index of regular variation equals —57/16. 

If the assumption E(X 4 7 ^ 8 e ) < oo is not satisfied, Theorem 3.3 may still be helpful. For instance, let us 
assume that X^ 1 is regularly varying with index —1/2, so that the above moment assumption does not hold. 
But since this assumption implies that A 4 ' is regularly varying with index —1, we can write the above 
joint exceedance probability as mU xf j > ax, 1 < i < 3) for Xj = Xj, 1 < j < 3, X A = X 4 1/2 and 

/ 1 -0.5 0 0 \ 

A = (ffl-,;j )l<i<3,l<j<4 = 0 1 —0.5 0 

\ 0 0 1 1 / 


If we replace A in (3.2) by A, then the optimal solution is given by k = (3/2,1,0,1) and this solution is 
unique and non-degenerate. Furthermore, 


(l J A^A) a 


/ 

1 T 

V 


1 -0.5 0 \ 1 \ 

0 1 Oj A = ((1,3/2,1)A) 3 = 1/4. 


Since E(X 3 +t ) < oo for e £ (0,3/4) by our above assumptions, we can apply Theorem 3.3 also in this case 
to obtain again that (Yi,Y 2 ,Y 3 ) is regularly varying on (0, oo) 3 with respect to scalar multiplication. But now 
the index of regular variation is equal to — ■_ 1 kj = —7/2 and the limit measure fi is given by 

A(x?=i(ci,oo)) = ^c/ 1 c^ 3/2 C3 1 S(X3 1/4 ). 
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4. Proofs and auxiliary results 

4-1. Proof of Theorem 3.1 

Proof. We start with the proof of (a). The optimal solution k to (3.2) lies in the closure of 
N{ A) := {z € R m : Az > 1} . 

Since the ray {z £ R m : z = (1 + 6)n,6 > 0} is a subset of the open set N( A), for all e > 0 there exists 
e' > 0 such that 


Qy + 

Thus, for x > 1, 




- + 1 


2Er=i«i 


+ e C N(A). 


P ^JJ Xj ij >x,l <i<n 

( m 

> P | Xj' J > x,l < i < n, and A ? >0,1 < j < rn 


U=i 


P (MM 

\ V ln(x) 


G iV(A) 


> 


ip 

i=i 


1 < m 

( !+oS'm . K . ) ~ e ' 


X 


9 r'm — I — t Kj ( 1+9 —— 

^3 = ! 3 / < Xj < X V 


+e' 


By the regular variation of Xi,..., X m , the expression on the right-hand side is of larger order than 

)+e'-. 


n 

3 =1 




= x -E”i«r|W-| > x -T.7 =1 nj-e t x > 1; 


which proves (a). 

For the proof of (b) let us for simplicity assume that c > 1 so that P(Xj > 1) = 1 for those 1 < j < m 
with Kj = 0. The modifications for general c > 0 (substitute Xj/c for Xj) are simple. Let A be as in Lemma 
4.2 (a) (see Section 4.3 below), so that we have 


J X°" 3 > x, 1 < i < n 

3 =1 
m 

=> Xf j > x, 1 <i<n 

3 = 1 
m 

=>■ ma x(Xj, l) aij > x, 1 < i < n, 

,?=i 

where we have used that > 0 for all 1 < j < m with Kj > 0 and Xj > 1 for 1 < j < m with Kj = 0. The 
last inequalities imply that for x > 1 


m 


y. + 

3=1 


ln(max(Xj, 1)) 
ln(x) 


> 1, 


1 < i < n 
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ln(max(Xj, 1)) ^ 

In (re) — ^' 


3 = 1 


"J ! 


J=1 


because otherwise k could not be an optimal solution to (4.36), in contrast to Lemma 4.2 (a) and our 
assumptions. Thus, for rr > 1, 


P ( Xj' 3 > x, l<*<nJ<Pj]^[ max(X 7 -. 1) > x^-o =1 Kj 


(4.1) 


yi=i 

By our assumptions, for all 8 £ (0,1), 

1-5N 


u'=i 


e 1 1 n rmix ( x j ,!) 

U'= 1 


= E (max(X 7 -, l) 1 s ) < 

3 =1 


oo 


and by the Markov inequality and (4.1) we conclude 
Pm X“ ij >x, 1 <i<n\ = 0 (x~ {j: T=i 


u=i 


for all (5 £ (0,1). Choosing 8 < e/ Y^jLi K j yields (3.5). 
4-2. Proof of Theorem 3.3 


□ 


In order to prove Theorem 3.3, we first deal with a setting that covers a slightly more general case for the 
solution of the linear program (3.2) than the one assumed in the statement of Theorem 3.3. The proof of 
this result is by induction on the number of positive components in the unique optimal solution to (3.2). 
Several auxiliary results needed for the proof can be found in Section 4.3. 

Proposition 4.1. Let A = (a^) £ R raxm w ffli n,m £ N be such that the solution k to the linear optimization 
problem (3.2) is unique with A k = 1. Define J = {j £ {1,..., m} : Kj > 0}. 

Let Xi,... ,X m be independent non-negative random variables. Assume that 

for j £ J : Xj is regularly varying with index —1, 
for j £ {1,..., to} \ J : P(Xj > 1) = 1 and E(Xj^ s ) < oo for all 8 £ (0,1). 

Then, 

p (n;ii x -' 3 > m < * < n) 

lim —-- 

rc—>-oo n PM > X K i) 

3&J 

= J II x j 2 M d (xj)jej) ® P {X}) ^ J Kxj)^) £ [0, oo) 

M( A) ■ jeJ 

with M( A) := {(aq,... ,x m ) £ [0, oo) m : ]lj=i ^ 1 1 <i <n}. 

Proof. The proof is by induction on the number l of positive components in the unique optimal solution 
k. Note that n > 0 and that at least one component of n has to be positive in order to satisfy A k > 1. 
In the following, we assume w.l.o.g. that the first l £ N components of n are positive and the last to — l 
components are equal to zero (if this is not the case, interchange the A}’s and the corresponding columns 
of A accordingly). 
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We start now with the case 1 = 1, i.e. k± > 0 and Kj = 0 for 2 < j < m. Our assumptions imply that 


djiKi = 1 for all 1 < i < n, i.e. an = ... = a n i and Ki = a a , 1 < i < n. Thus, 


P X T‘ >*• 1 < i < » 1 = p U, mm (jj X. 
If the linear program 


Q>ij ^1 

3 


> x K 


find x £ [0, oo) m 1 such that 


has no feasible solution, then 


a 12 


min y aijXj-i <1, Vxg [0,oo) 

l<i<n ^' 


m —1 


3=2 


and thus 


an 


\ a n 2 ■ * * a,. 


X > 1 , E Xi — > min! 


(4.2) 


(4.3) 


mm 

1 <i<n 


U=2 


< e 1 a.s., 


(4.4) 


because ln(X 2 ),..., ln(X m ) > 0 almost surely by our assumptions. On the other hand, if there exists 
a feasible solution to (4.3), then there exists e > 0 such that all feasible solutions x to (4.3) satisfy 


JO™ / x i > K l + since otherwise there would exist a solution x' = (0, x[,..., x' m _i) T («i, 0,..., 0) T to 
(3.2) with JO = i a/- < in contradiction to our assumptions. Hence, an optimal solution to (4.3) exists 


with JO"0/ Xi > Ki + e. By Theorem 3.1 (b) we have 

n*H >x ) = o (a; _1 " e/(2Kl) ) 


P | min 

l<i<n 


X —> 00. 


U= 2 


So, whether there exists a solution to (4.3) or not, we have 


E 


( 


\ 


1+€/(4«i)N 


mm 

l<i<n 


n* 

\ 3 = 2 


aij k i 

3 


< OO 


and we may thus apply Breiman’s Lemma, cf. Breiman (1965), to derive the asymptotic behavior of (4.2) 
as x —oo. This gives us 


lim ■ 

X — YCC) 


p (n ?=1 X T* > *, 1 < < < «) 


P(X i > x Kl ) 


= E min TT X 

\ 1 <i<n I 11 

\i = 2 


^1 
3 


x x 2 l(da;i) <g> p( x ih<j<m (d(xj) 2 <j<m) 


{x6[0,oo) m :ji>maxi< i <„ (iIJLa x j } 
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/, 

{xe[ 0 ,oo) m :njLi x j J >l,l<i<n} 


X 1 2 l(dxi) <g> p( X ih<i<™ (d(Xj)2<j<m), 


which concludes the proof in the case 1 = 1. 

For the induction step, assume that Proposition 4.1 holds for all matrices A* G K" xm ,n*,m* G N, for 
which the corresponding linear program (3.2) (with A replaced by A*) has a unique solution k* and for 
which A*k* = 1 and at most l — 1 > 1 components of n* are positive. In the following, assume that n has 
l positive components, again w.l.o.g. the first l ones. 

Define the map <8> K as in Section 2. From Lemma 2.5 we get that p( A b) i<i<™ j s regularly varying on 0„ 
with respect to ® K . 

Now, with An = 1 we have 


P 


(ft X T j >x, 1 < i <n^j 


P 



> 1 , 1 < i < n 


P(x M), 


where X = (Ar, ..., X m ) and 

{ m 

x G [0, oo) m : n X a / j >1, 1 < i < n 
l=i 


For S > 0 write 


P(X £x® k M) _ P (X G x (M n ((J, oo) 1 x [0, oo) m - 1 ))) 

nU p(Xi > ^" nUi p (x > xKi ) 

P (x G X ® K (M n {(6, oo) 1 x [0, oo) m - l ) c )) 

YlLiP( x i> xKi ) 


(4.5) 


We will first show that the second summand in (4.5) tends to zero as first x —> oo and then S \ 0. Note 
that 


< 


lim lim sup 

<5\0 x— xx) 


P(Xgi® k (Mn ((5, ooy x [0, oo) m - ; ) c )) 

nU p&i > 


l 

) lim lim sup 

s\0 x^oo 


P(X G x (g) K M, X k < 5x Kk ) 

I?i = iP(Xi>x*<) 


(4.6) 


We will show that all summands in (4.6) equal zero. To this end, note first that we may apply Lemma 4.2 
(b) to the matrix A, i.e. there exists a matrix A such that k as above is the unique solution to the linear 
program (4.36) with aij > 0 for 1 < i < n and 1 < j < l and Ak = 1. We have 


P (X G x M, X k < Sx Kk ) = 


< 


P 


^fl Xj iJ >x, 1 < * < n, X k < Sx Kk j 


P ( Xf j >x, 1 <i<n, X k < Sx Kk 
U=i 


(4.7) 
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by Lemma 4.2 (b). For ease of notation, we restrict ourselves to the analysis for the summand k = 1 in (4.6). 
For C > 0, use (4.7), an > 0,1 < i < n, and A/t = 1 to write 

P(X e M,X 1 < 6x K ') 


< 


< 


(n”i K 


P (n?=i X T 3 > x > l<i<n,X x < 


nu p{Xi> x *<) 


I P{£ r >z~\X 1 < 

[0,oo) 


m in ff (J^y ij/aii 

1 <Sx Kl )P 1 ^ n ^ Kx 3J (dz) 


nUi p(x l>x ^) 


min T 

/ P(-§k> -z” 1 ) / p ( x 1 > x' 41 )^ 1 - 4 -"^ 2 

(,5- 1 ,a:''i/C] 

n l i= 2 p(x i >x^) 

min J 

/ P(£ r >z~ 1 ,X 1 <6x K ')P 1 ^ n t= 2 

(x K i /C, oo) 


x i \ a ij/ a H 


X i \ a ij/ a i. 


(dz) 


(dz) 


+ 


nu p&i > *“* 


=: I(x,S,C) + II(x,S,C). 


(4.8) 


We deal hrst with I = I(x, S , C). Use aij/ctn > 0,1 < i < n, 2 < j < l, and k 3 = 0 and P(Xj > 1) = 1, l < 
j < m, to obtain 


mm 

l<i<n z ' a,l v x 

- - j=2 V 


in y in ( *1) < min V 

<r * J n 1 V r K n I '\<T'i<''n * ^ 


l<i<n L ^ QjA i 
- - j=2 


In 




a.s. 


(4.9) 


Choose e £ (0,1) according to Lemma 4.3 such that the expression on the right hand side of (4.9) is a.s. 
bounded by 


1 — e 

1+7 


1=2 


E( ln (^j 


Xi 


For this e > 0, there exists C > 0 such that 


P(Xi > x Kl ) 


<(l + e)z 1+e VI < z < x Kl /C, x > C 


by Potter’s bounds applied to x ^ P(Xi > x) (cf. Bingham et al. (1987), Theorem 1.5.6). So, for S < 1, 
the numerator of I(x, S , C) is bounded by 


’ (<5 1 ,oo) 


min n ( 

(1 + e )z 1 + e P 1 - i - n 3 = 2'' x 3 ' ( d ~) 


M( A,5) 


(1 + e) min TT z\ 

l<i<n J-J- ■> 


(l+e)ay/a«i 


1=2 


< 


(1 X e ) IT max(l, Zj) 1 ~ e P' 

M(A,5) j_ 2 


J 2 <^<m (dz) 
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< 


tfih 


E / 

^^{0,1—e} m_1 J M(A.,8) 


(i + e) n 


3 =2 


x 1 ^2<3< m (dz) 


with 


M( A, 6 ) := <( (z 2 , ■ • •, z m ) G [0, oo) m 1 : min TT z ^ /atl > $ 

l<i<n J 

~ ~ 3=2 

Note that y^”' ,, dik^k > 0,1 < i < n, by our assumptions about A and k and let 


D(5) := min 5 a ^l^k=i a ^k > g 

l<z<n 


(4.10) 


and 




<2<n,2<j<m — 


E m ~ 

fc— 2 Q'ikK'k 


l<z<n,2<7<ra 


Hence, 


M(A,<5) C 


m 

(^ 2 , • • •, 2m) € [0, oo) m_1 : min TT z“* J > £>((5) > 

l<z<n A A J 
-- 3=2 J 


(4.11) 


Note that a feasible solution to the linear program 

m 

find x = (&2 , ■ • ■ ,x m ) T > 0 such that A'x > 1, -A- min! (4-12) 

i=2 


is given by = (« 2 , ■ ■ ■, n m ) T with A'n' = 1. Furthermore, this is also the unique optimal solution to 
(4.12), because if there would be another feasible solution (x 2 , ■ ■ ■ , x m ) T to it with Y^j= 2 x 3 — SJI 2 
then x' := (ki,x 2 , ... ,x m ) T would be a solution to (4.36) as well because of 


Ax' = 


Oii«i 


E«-' 

3=2 


Xj 


> 


ClaK 1 


1 <i<r 



= 1, 


Ki<n / 


as A/t = 1. This would lead to a contradiction to our assumption about the uniqueness of k and Lemma 
4.2 (b). Thus, 


' ^ E 

(/£ )2<j<md0il e} m 1 


E 

(Pi)2<j<m£{0,l— e} m 1 


z: min n zA ] >D{5) 


l<i<n / * 


(i+e) n ;= 2 dz) 


num.>^) 


/ 


z: min n >1 


l<i<n 


nr= 2 <^ 


\(D(S)*1 3 


3<m(dz) 


nL^> w)*) re< ) 
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n,.2 p(*i > 


j =2 


Now, by Proposition 4.4 combined with the induction hypothesis and the properties of A', the first factor 
of each summand in the above expression converges to the finite expression 


/ m l 

n n x 7 2 m*jh<j<i) ® 

M(A') j=2 j=2 

as x —»• oo, whereas the remainder of the expression converges to 
i 

(l + e)II W' s(1_/3j) 

i=2 

with Kj (1 — ) > 0 for 2 < j < l by our assumptions. The first limit does not depend on the value of <5 > 0, 

while the second converges to 0 as <5 \ 0 and thus D(S) —> oo by (4.10). We have thus shown that 

lim lim sup I(x, 6, C) = 0 

x —^oo 

for C large enough. Let us now deal with II = II(x 1 6 , C ) from (4.8). We have 

. aij j flj 


I '**' / V \ Q'ij / 


II < 


nL p(x i>x *<) 


< 


< 


where we set 


P J] max(l ,Xj)^ > x/C 5il , 1 < i < n 

v =2 _ l 

ni=i -P(max(l,Xj) > x Ki ) 

( m 

P ( II max(l, Xj)“y > xD'(C), 1 < i < n 

<3 =2 


IL=i P(max(l, Xi) > x Ki ) 


n! = i -P(max(l, Xj) > x Ki ) 
n l i=1 P(X i >x*‘) 

> x Ki ) 


(4.13) 


D'(C) := min C~ hil > 0 

l<i<n 

for abbreviation. Set 

A" := {a,ij)\<i< n 2<j<m € ]R nx ( m_1 ) 
and consider the linear program 


m 

find x = (x 2 , • • •, x m ) T > 0 such that A^x > 1, Xj —> min! (4-14) 

i—2 


If this linear program has no feasible solution then 


m 

min IT max(l,X,-)° w < e a.s. 

KiKn 1 - 1 - 
~ - j 2 
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(cf. (4.4) for analogous reasoning), so the first factor in (4.13) equals 0 for x large enough. On the other 
hand, if there exists a feasible solution to (4.14), then there exists e > 0 such that all feasible solutions 
(x 2 ,...,x m ) T to (4.14) satisfy Y^j= 2 x o > J2jLi K j + e, since otherwise there would exist a solution x' 
= (0, X 21 ■ • •, x m ) T ^ k to (4.36) with (Cjli x j — Xqli K j; in contradiction to our assumptions and Lemma 
4.2 (b). In the latter case, the numerator of the first factor in (4.13) is of smaller order than x k 3 _£ / 2 
as x —> oo by Theorem 3.1 (b), while the denominator is regularly varying in x with index — Xq=i K j an< ^ 
the second factor in (4.13) equals 1 for x > 1. So, in both cases, and for all C > 0 


lim limsup II{x, (5, C) = limsup II(x, <5, C) = 0, 
s\0 x—>oo x—±oo 


and the first summand in (4.6) is equal to zero. All other summands can be treated analogously. 
Taken together, we have shown that 


lim lim sup 

x —>oo 


P (X £ x ® K (M fl ((<5, oo) 1 x [0, oo) m l ) L )) 

nUl TO > *"*) 


= 0. 


With //(•) as defined in (2.2) we have 


(4.15) 


P 


lim 


(nr=i X T 3 > x ^<i<n) P(X £i» K M) 

= lim 


nL p i x i > xKi ) s ^° c nu p ( x i > xKi 


= lim lim 

x^oo 


F(Xei® K (Mn ((5, oo) 1 x [0,oo)"*-'))) 

nUi to > * Ki ) 


= lim p (M fl ((5, oo) 1 x [0, oo) m *))) 


= lim 

S\0 


Mn((6,oo) 1 x[0,oo ) m-1 ) 1 


x t 2 \(d(xj)i<j<i) ® p( x Pi<i<m(d(xj)i < j< m ) 

3 =i 

J n X T 2 Md( x j)l<j<l) ® (d(Xj-)/<j<m), 


3=1 


by monotone convergence where we used that 
M n ((0, oo) 1 x [0, oo) m ~ l ) = M, 

because Kj > 0 for \< j < l implies that at least one 1 < i < n exists with > 0. But then 

mini<j< n n™i X^ il = 0 as soon as Xj = 0,1 < j < l, so M C ((0, oo)* x [0, oo) m_i ). 

This concludes the proof of Proposition 4.1. □ 


Proof of Theorem 3.3. Let us again assume w.l.o.g. that the first n > 1 components of k are positive and 
the last m — n> 0 components are equal to zero. We start with some implications of our assumptions about 
the matrix A. Since the optimal solution k to (3.2) is unique it must be a vertex of the polygon defined by 
{x > 0 : Ax > 1}, cf. Sierksma (1996), Theorem 1.5. Each vertex of {x > 0 : Ax > 1} corresponds to a 
so-called basic feasible solution (cf. Sierksma (1996), Theorem 1.2) of the standard form linear program 

m 

find x £ R m+ra such that (A; (—1) ■ E„)x =1, x > 0, Xi —> min!, (4.16) 

i=i 

where the matrix (A, (—1) • E„) £ H nx ( m + n ) consists of the columns of A in its first m columns and of the 
columns of the n-dimensional unit matrix, E n , multiplied with —1 in its last n columns. The basic feasible 


imsart ver. 2014/10/16 file: randomproductsjanssendreesarxivv2.tex date: May 13, 2016 






Joint exceedances of random products 


17 


solutions of (4.16) can be found by choosing n linearly independent columns of (A; (—1) ■ E„) with indices 
B C {1,..., m + n}, denoting the resulting matrix by (A; (—1) ■ E„)b and deriving s b = ((sj)jeB) T := 
((A; (—1) • E„)s) _1 l. If sb > 0, then we call 


xs = (xi,..., x m+n ) T with 


Xj = 0 , 


if j e B, 

if j € (1,... ,m + n} \ B 


a basic feasible solution to (4.16). The corresponding solution to (3.2) is given by the first m components 
of xb, the remaining last n components of xg are called slack variables. Since we assumed that the first n 
components of k are positive and that the optimal solution is unique, it can only correspond to the basic 
feasible solution with B = (1,..., n} which implies that A K , the matrix which consists of only the first n 
columns of A, is invertible and (ki, ..., n n ) T = (A K ) -1 1 which leads to Ak = 1. Thus, the assumptions 
about A of Theorem 3.3 are a special case of the assumptions about A of Proposition 4.1. Furthermore, the 
optimal value of (3.2) equals 


2."i = 2^i = 1 ( A « )i- 

i=i 1=1 


(4.17) 


Since we assumed k to be unique and non-degenerate, the optimal solution to the dual problem (3.9) is 
unique and non-degenerate as well, cf. Sierksma (1996), Theorem 2.11. Furthermore, the optimal solution k 
to (3.9) is in our case given by 


k={A~ 1 ) T l, 


(4.18) 


cf. Sierksma (1996), Theorem 2.2. This explains Remark 3.5 (a). 

Again w.l.o.g. assume in the following that (l 7 A f 7 1 A)j = 0 for n < j < n' with n < n' < m and that 
(l 7 A~ 1 A) J - ^ 0 for n' < j < m. Define now for 1 < j < m 


Xj, 

for j < n , 


x h T _ 

for n < j < nl 1 


v( lT A re 1 A) i +e 

X j > 

for j > n', (l T A“ 1 A) i 

V 

o 

Y (l T A~ 3 A) j -e 

X J 

for j > n', (l r A- 1 A) J - 

A 

o 


Xj = 


with e > 0 as in the statement of Theorem 3.3 and set furthermore A with 


(4.19) 


an = 


MO , 


(1 T A^ 1 A)j — e ’ 

Obviously, this leads to 


for 1 < i < n,j < n, 
for 1 < i < n,n < j < n', 

(i T A^A) J+ e ’ for 1 < i < n, j > n', (l T A~ 1 A) j > 0, 
for 1 < i < n, j > n', (l 7 A k : 1 A) j - < 0. 


(4.20) 



> CiX, 1 < i < n \ = P 



> CiX , 1 < i < n 


(4.21) 


In order to apply Proposition 4.1 we will distinguish between events where Xj > 1 and those where Xj < 1 
for n < j < n'. Therefore, for J C {n + 1,..., n'} denote the event {Xj > 1 for j £ {n + 1,..., n'} \ J, Xj < 
1 for j £ J} by B(J), and write 


P 


Xj' 3 > dx , 1 < i < 
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^ P [ ]^[ Xp 3 > ax, 1 < i < n 


JC{n+l,...,n'} \j=l 


B(J) P(B(J)). 


For J with P(B(J )) > 0 define now independent random variables Xp\ 1 < j < to, with 


(4.22) 


P x< i J) = <j p x i 1 l^< 1 j for j € j ; 

for j G {n + 1, ..., n'} \ J. 

Furthermore, set with 


' P Xj , for j G n.n' + l,..., to}, 


-(j) = I a*j, for 1 < i < n,j G {1 

I "®ij) for 1 < i < n,j G J. 


(4.23) 


(4.24) 


By independence of the Ay’s and of the the first factor of each summand in (4.22) is equal to 


m - (J) 

n (*1T 


> CiX, 1 < i < n 


(4.25) 


u=i 


for all J C {n + 1,..., n'} with P(B(J )) > 0. The vector n is a basic feasible solution to 

m 

find x > 0 such that A ( J ^x > 1, —> min!, 

i=i 

because A IJ> k = 1, as the first n columns of A ( ' J ' > are identical to those of A. Furthermore, 

f(l T A-'A)j = l, 
e _1 (l T A~ 1 A) 7 - = 0, 
e-\l T A-\-A))j=C 
0 , 1 ), 


(4.26) 


if j 

VI 

n. 




if j 

G 

{«+!,••• 


\J 

if j 

G 

J 




if j 

> 

n ' 

,(1 T A k 

'A), 

o 

A 

if j 

> 

n / 

,(1 T A k 

'AJj 

A 

o 


(4.27) 


which proves that k, is the unique optimal solution to (4.26), because 1 — (l 7 A^Af 7 ))^ is strictly positive 
for all non-basic variables n < j < m. (cf. the analogue of Theorem 1.6 and the remark after the proof of 
this theorem in Sierksma (1996) for a linear minimization problem instead of a maximization problem). 

Let now J C {n + 1,... , n'} and let A l - J ' ) be the matrix described in Lemma 4.2 (c), corresponding to 
with a[’p > 0 for all 1 < * < n, j G {1,..., n, n’ + 1,..., m}. Then, for c > 0, 


=P 


+ P 


m 

a (J) 

\ 

n W) 

a ij 

> ax, 1 < * < n 1 

0 = 1 / 


a (J) 


n (*f ,> ) 

ij 

> ax, 1 < i < n, Xj J ' > 

0-1 

' m 



ii « 

P 

> ax, 1 < i < n, 3 n! 


pj) 


(4.28) 


U=i 
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and for the second summand we have 


P J > ax, 1 < i < n, 3 n' < j < m : X^ < c 


U=i 


z<.j) 


< 


y Wn(*JT> 

$^Kc{n' + \j =1 


x mm Cfc, 

l<k<n 


x f ] <c,j G K, X] J) > c,j e {n' + 1,... ,m} \ K 


3 

~(J) 


< 


£ W n (*1T n 

0 ^K'C{n / + l,...,m} \ j£{l,...,n'} j^{n' 




. ~ £ a 

> x min Cfc min c 

l<k<n \ l<i<n 


(■/) N 


Set 


1 < i < n . 


^ -( j )' 

— > a. 


min c > 0 

l<i<n / 


D(c) = D(c, ci,...,c n ) := ( min c fc I 

\l<fc<n / 

for abbreviation and note that by our assumptions and Lemma 4.2 (c) k is the unique optimal solution to 

~ 771 

find x > 0 such that A (,77 x > 1, Xj 


lint 


i= 1 


with A^k = 1, that P x t = P x t , 1 < j < n, that X^ x ,..., X^p > 1 a.s. with E(X^) < oo, n < j < n', 

and that E( max(l, xj' 7 ' 1 )) < oo,n' < j < to. Therefore, apply Proposition 4.1 to obtain for 0 ^ K 
{n' + 1,..., to} 

;(J) x(J) 

p( II (xP') %3 II max (l. t3 > D(c)x, 1 < i < n) 

V je{l,...,n'} V 7 je{n'+l,...,m}\K V 7 7 


c 


rii=i p ( x 3 > x Kj ) 


-)• (D(c))" e ?=i K3 D(J), 


as x —> oo for some finite constant D( J) which does not depend on c. As D(c ) —> oo for c 0 we conclude 
from (4.28) that 


P 


lim 


m j) f 


> ax, 1 < i < n 


n =i p{Xj > x-3) 


-M) 


P 


= lim lim ■ 

c\jO x—>oo 


J} (xj J ^ 3 > ax, 1 < i < n, Xj' 7 ' 1 > c,vl < j < to 

: nu p(*3 > 


(4.29) 


Let now 


v - 'm - (■/) v—'m 

Cj(c) = ac~ “o' = CjC~ 2 -'i=»'+i “ y , 1 < i < n. 
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Then 


P 


lim ■ 

x—>oo 


n iJ > CiX , 1 < i < n, Xj' 7 ^ > c,n' < j < m 

n;u p(Xj > x*<) 


= n nx\ j) >c) n^>^) 


j=rt' + l 


Vj=i 


/ 


~(J) m 


n(*r)“ ii 


x 


{j)\ °« 


> Ci(c)x, 1 < i < n 


ij=1 2 =n'+l \ / 

The last factor in the above expression can be written as 

~(J) 

{3=1 

where xj J,c \ 1 < j < to, denote independent random variables with 

x (J) 


min > c 


p 


m - ( J) 

]-[ ^xj J ’ c )j ' 3 > Ci{c)x, 1 < i < 


(J.O) P A o 


P x i ' = 


for 1 < j < n', 

pc X j X j > C) f Qr n ' < j < TO. 


Set 




nilic^c)^ 1 )-, 1 <j<n, 

n < j < to, 


which implies that 


1 Cj(c) a « = exp I ^ ay ^](A K 1 ) jfc ln(c fe (c)) I = Ci (c), 1 < i < 


n n 


3= 1 

Then we have 

P 


3=1 


ki=i l 


m / , T alT 

n( j r) 


> Ci(c)ir, 1 < z < 


nL^>^) 


n n 


Cj (c) 




m / x' ij 


> x, 1 < i < n 


n"=iP(^>^(c)*"o 


n-=i > Cj(c)ir^) n'~i -P(^J > ) 

As a: — > oo, in view of (4.18), the second factor of (4.32) converges to 

n n n n n 

n&to )- 1 =n ii(c*(c))- (A --^ ==n 

j— 1 j = 1 2=1 2=1 2=1 


(4.30) 


(4.31) 


(4.32) 
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n 

n( 


C 7 k * C ki + 1 a ij’ ) = C E"= i £" ln'+l S « ) « 


n< 


Note that P^ = P Xj for 1 < j < n and 


X^ c ) / ^( J ’ c ) \ 

—--= > 1 a.s., E [ J ] < oo for n < j < m. 


fo'(c) 


fo(c) 


Hence, we can apply Proposition 4.1 to see that the first factor of (4.32) converges to 

r. n 

/ TT xj 2 \(d(x j ) 1 <j<n) <8> P ( ^ J< ’ (d(Xj) n< j<m) 

Jm( A«)“ 

Write xi = (aq,..., x n ), X 2 = (x n +i,..., x m ) for abbreviation and use the substition y = 
(ln(xi),..., ln(x„)) to see that the above expression equals 


[OjOo) 71 


/ n 

ex P( _ E J/i )*( dy )- p(Xi< ’ )) " <3 '- m (dx 2 ) 
1 i=1 


yeH"-:exp(A K y)>(njl n+ i x j ' 3 )i <i< r , 


[0,oo) 


/ 


AJ) 


det(A K )| 1 

6R":exp(z)>(n”L„ + i *} '* )l<«<n| 

n n \ 

exp ( - )kizi J l(dz)P (A i ’ ) )™<i<m(dx 2 ) 


\ k= 1 1 = 1 

= | det(A K )|~ 1 J 

[0,oo) m_ 

| det(A K )| _1 

T n 


n 

1=1 


exp ( -*,£( Aj) kl ) &ziP ( 'X?' C) ' >n<i ^ m (dx 2 ) 


a ( J ) 


k = 1 


(ln(nr=n+i x j 3 ),°°) 

n m ^ (J) 


rii=l((^K ) T l)i d[0,oo) m -” ;=lj - =rl _ | _ 1 

|det(A K )|- 1 ^ 

nr=i « 


n ri A ii ^ i(a ‘' 1 )w ^ m )»«»( dx 2 ) 


n £ (w) 


Er=ia i ( / ) ((A- i ) T i) i 


j=n+l 


I det(A K )| 

n?=i«i 


-1 


II H (W 1 ) 


v—»n ~ (JO ~ 


j=n'+l 


E (J) > c 1 c 


-Er=i a if«< 


= ( iTa « 1j 

(4.33) and (4.34) to see that the expression in (4.32) converges to 


where we used in the final step that = (1 T A K 1 A( J ))i = 0 for n < j < n 1 , cf. (4.27) 


idet(A K )i- i nr=iE i m 


n_i*< 


n e (*; j) ) 


EIU 


- ( J) , 


E 7) >c 


j=n '+1 

as x — 00 . Now, (4.35) together with (4.29) and (4.30) yields that 

p(n(.if>r 

lim —— 

x—>oo 


> CiX, 1 < i < n 


n;=i 


(4.33) 


ln(xi) := 


(4.34) 
. Combine 

(4.35) 


imsart ver. 2014/10/16 file: randomproductsjanssendreesarxivv2.tex date: May 13, 2016 











22 


A. Janfien and H. Drees 


= lim 

c\Q 


i detcA ^)!- 1 nr=i ■ 


— 1 TT n — Ki 171 


n * c 


x. 


( j )\ 


£?= i 


~(j) ■ 


rr= 


i=1 


j=n '+1 


m / 

n * c 


x. 


r) 




nr =1 ^ 

det(A K )i- i nr=icr: m 

j=n' + l 
m 

[| £?(xp= i04i **), 




idet(A K )|- i nr=icr: m 

i=™'+i 


rr= 


2=1 ^ 


where we used Ei=i®ij' > ** = ((A^) T (A“ 1 )’ 1 l)j = (l 7 A~ 1 A^ J ^)j > 0,n' < j < m (cf. (4.27)) in the 
penultimate equality and (4.19), (4.20), (4.23) and (4.24) in the final equality. This expression no longer 
depends on J C {n + 1,... n 1 } and therefore (4.21), (4.22) and (4.25) lead to 


P 


lim 

x—>oo 


(nr.. y u 


> CiX , 1 < i < n 


) 


P 


n"=i p(Xj > ) 

det(A K )i- i nr=icr 


= lim 


(nr.Ar 


> CiX , 1 < i < n) 


uu p ( x j > xKj ) 


nr=i«i 


= | det A k | 


, 11 ". 


j=n'+l 

-(i t a- 1 ) 


E(x'p*= iaiiki } 


2=1 °2 


n’UlUA; 1 , 


k- n E(y iTA; ' A) '). 


j:itj=0 


and so the limit in (3.7) equals the expression in (3.8). By Theorems 2.4 and 2.5 of Lindskog et al. (2014), 
this shows that c{x)P ^(aP 1 TIj=i ^j' J )i<i<n € , x > 0, with c(x ) = (Tlj.i p ( x j > x Kj ))~ 1 , is relatively 

compact in M(o i00 )n. Furthermore, all accumulation points of this family agree on a generating 7r-system. 

Thus, P iU T= i V)i<‘<™ is regularly varying on (0,oo) n w.r.t. scalar multiplication, cf. Example 2.4. The 
index of regular variation follows from Lemma and Definition 2.2 since c is regularly varying with index 
- EU = - E7= 1 = -! TA cf. (4.17). ' ' □ 

^.5. Auxiliary results 


In the following, we collect two lemmas and a proposition which are needed for the proofs in Sections 4.1 
and 4.2. 

Lemma 4.2. Let n = (sq,..., n m ) T be an optimal solution to (3.2). 

(a) There exists a matrix A = (a^) € R nxm such that 

— the columns j in A for which > 0 have all positive entries, 

— k is an optimal solution to the linear program 

m 

find x > 0 such that Ax > 1, 'y^ Xj —> min! (4.36) 

2=1 

for all x, Xi,X 2 , ■ • ■, x m > 0, 

m m 

x a P > x, 1 < * < n => xf j >x, 1 < i < n. (4.37) 

3 =1 3=1 

(b) Moreover, if the assumptions of Proposition 4-1 hold, then the matrix A can be chosen such that 
additionally 
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— k is the unique optimal solution to the linear program (4.36), 

— An = 1. 

(c) If the assumptions of Theorem 3.3 hold, then there exists a matrix A = ( dij ) £ K" xm SMC /j 

— the columns j in A for which (1 T A~ 1 A)j > 0 have all positive entries, 

— k is the unique optimal solution to the linear program (4.36), 

An = 1, 

— for all x, Xi,X 2 , ■ • ■, x m > 0 (4.37) holds. 


Proof. First note that if > 0 for all 1 < i < n and all j such that Kj > 0 (cases (a) and (b)) or 
(l T A^ 1 A)j > 0 (case (c)), then we may simply set A = A. So, assume the contrary in the following. Set 
J := {j £ {1,... ,m} : Kj > 0}. Since we have assumed an optimal solution k to (3.2), there also exists 
an optimal (not necessarily unique) solution k = (ki,..., k n ) T to the dual problem (3.9) and this solution 
satisfies ^"=i = Xqli K :i- cf- Theorem 2.2 in Sierksma (1996). Furthermore, by the Complementary 

Slackness Theorem (cf. Sierksma (1996), Theorem 2.4) we have (A T k)j = 1 for all j £ J. For assertions (a) 
and (b) let a m i n := —(mini <i<n,jej otij) + e for some e > 0. By our assumptions, a m i n is positive. Define 


A — (Aij )l<i<n,l<j<m with 


j 1 Amin 2-^k— 1 Q'kj^'k 
1 i A m iri 2^k—l K'k 


As seen above, we have 'fZ'f.-i a kjkk = 1 and thus dij > 0 for j £ J and all 1 < i < n. 
Note that 


A.K — ^1 -{- Qmin ttij ( {o.ij ~\~ ttmin) l<i<n,j£ j) ((^bOy'Gj) 

/ m \ — ^ / m \ 

^ I 1 T U m i n ^ ) Hi J I 1 T U m i n 'y ' re^l I = 1, 


(4.38) 


so k is a feasible solution to (4.36). Furthermore, if there would exist a n' > 0 with A k,' > 1 and K i < 

YJiLl then 


E ( T a min E akjkk J Kj- > 1 + Amin E Kfc, !<*<«, 

j —1 V k— 1 / k— 1 


and thus 


E a b K 7 - 1 + Amin ^ ^ ^k A m i n EE Ofcj ii k 


i=i 


fe=l 


fc=l j = 1 




m m 


^ 1 T «min I ^ ^ 

\fc=l .7=1 

> 1, 1 < j < n, 


(4.39) 


(4.40) 


where we used in the penultimate inequality that X!fe=i a kjkk < 1 and re' > 0,1 < j < m. But this implies 
that k/ with K i < would also be a feasible solution to (3.2), in contrast to the assumption 

about the optimality of h. Thus, h is also an optimal solution to (4.36). 
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We are thus left to show (4.37) for the proof of (a). For x, x \, X 2 , ■ ■ ■, x m > 0 such that ]~f '=1 x^ 11 > x,l < 
i < n, we have 


n 


a ij a 

X 3 


min 


i=i 


rj^a m i n K,i 


1 < j < n, 


with strict inequality if ki > 0, which must be the case for at least one 1 < * < n. So by multiplication of 
left hand sides and right hand sides we obtain 






O- 


~[ X T n II x] kjaminkk > x x a ^nfi k = X 1 +a ml „££ =1 


u =1 


k =1 i =1 


fc=l 


J'J' a , a y+ ami n s, jl+Omin Efc=l 1 < * < 


J = 1 


TTl S~ a min Sfc—i a kj^k 

X ' 1+ “”in E ^ 1 ** > x, 1 < i < n 


j=i 


^ x^' 3 > x, 1 <i<n. 
i =i 


1 < i < n 


Thus, (4.37) holds. 

For the proof of (b), we use that the additional assumption implies that k is the unique optimal solution 
to (3.2) and that A/t = 1. Similar to (4.38) one shows that Ak = 1. Furthermore, if there would exist a 
k' ^ n with Ak' > 1 and Y^'p=\ K 'j — l K j then one shows analogously to (4.40) that this would imply 
that the optimal solution k to (3.2) is not unique. This shows that k, is the unique optimal solution to (4.36) 
and proves (b). 

For the proof of (c), we use that the additional assumption implies that k is the unique optimal solution 
to (3.2) and that k. = (A” 1 ) 7 ! is the unique solution to (3.9), cf. the beginning of the proof of Theorem 
3.3. Let for some e > 0 


a 


(c) 

min 


nun -f-j- b e 

l<i<n,i:(l r A- 1 A) ;l >0 {l 1 A k A)j 

o^_ +£ 
— i &kj ftk 




i V '' 71 

A-k. A)j >0 Z_^k 


which is positive by our assumptions. Define 


A (c) = (oj?) 


l<z<n,l<j<m 


with da = 


_ I „(c) V— 

a ij + a m j n 2^k=1 a kjKk 

1 I „( c ) V~' m 
1 + a min zJfc=l K k 


We have thus > 0 for those j with (l 7 A K 1 A) J > 0 and all 1 < i < n. The rest of the proof for assertion 
(c) follows analogously to the proof of (a) and (b) which did not depend on the value of a m i n >0. □ 


Lemma 4.3. Let the assumptions of Proposition J^.l hold and assume in addition that aij > 0 for all 
1 < i < n and those 1 < j < m for which Kj >0. Then for all j with Kj > 0 there exists e > 0 such that 

min V' — x k < (1 — e) V' x k (4-41) 

l<i<n did 

1 <k<m,k^j 

for all (xfc)i <k< m ,kjtj G [0, oo) m_1 . 
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Proof. For ease of notation and w.l.o.g., let us assume that Ki > 0 and treat only the case j = 1. For 
{xk) 2 <k<m = 0 the inequality holds for all e > 0. The rest of the proof is by contradiction. Let (e;);eN 
be a sequence such that ei > 0 for all l £ N and ei \ 0 for l —> oo. Assume that for each l there exists 
( x < 'f) 2 <k<m € [0, oo) m_1 \ { 0 } such that 


m 


mm ~ x \ > (1 - e;) 
\<i<n^ at 1 3 ^ 

3 =2 


„(/) 


k= 2 




m 


X< 

3 = 2 


J) 


‘y (0 - ^ 


1 < ? < n. 


E m 

k— 


fc=2 •''fc 


where we used that Oji > 0,1 < i < n, by our assumption. Define now 




„(0 


* ( ° = (^ ) ,...,kW) t = (o,«2+ * «i, •••,«>»+ (0 '«i I >0 

V E fc =2 4 E fe = 2 *fc / 

for all 1 g N. We have 

m m 

j=i j=i 

for all l £ N. Furthermore, 


(4.42) 


^ ' da n 


~{l) 

*3*3 


— X] a *J ( + 


.(0 


i=i 


J=2 \ 

m m 

= yy dij Kj + 


~(0 
Z^fe=2 x k 


«i 


i=2 


xf 

W) Kl 

j=2 Z-//c=2 x /c 


> 1 — ajiKi + (1 - ez)aiiKi = 1 - ejOiiKi, 


(4.43) 


for all l £ N, where we used Ak > 1 and (4.42) in the last step. For l —> oo, the bounded sequence must 

have an accumulation point k ^ k (because iti = 0 < «i) and /« > 0. But 

mm m 

yy e= x ** and x ^*3 > 1 

j=i i=i i=i 

by (4.43), so our optimal solution k would not be unique, in contradiction to our assumptions. Thus, for 
some e > 0, the inequality (4.41) holds for all ( Xj) 2 <j< m G [0,oo) m_1 . □ 


Proposition 4.4. Assume that 

lim P (PQ)i<j< m e y ® K M( A)) 

y->°° II P{Xj > y K i) 

j = 0 

= J n xJ 2 Md(xj) {j:Kj>0} ) ® P^^J^ (d(xj) {j:K}=0} ) G [0,oo) (4.44) 

M(A) ^ :k 3 >0 

holds for all X\, ..., A m and all matrices A which satisfy the assumptions of Proposition 4-1, with k being 
the unique solution to (3.2) and M( A) as in Proposition 4-1. Then also 

I nr=i x? i P( x ih Ki hi<i<'»} (dx) 

M(A) 

lim - 

y—^oo n P(X j >y">) 

j:Kj>0 
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/ m 

II II X j 2 H d ( x j){j.nj>0}) ® /’ 01 (d(Xj) {j:Kj = 0}) 6 [0, oo) 

.7=1 


(4.45) 


M(A) 


/or all Xi ,..., X m , A and k as above and all (3j G [0,1), 1 < j < m. 

Remark 4.5. As the preceding proposition is used in the induction step of the proof of Proposition 4.1, 
the convergence (4.44) had to be assumed. However, since Proposition 4.1 shows that (4.44) holds for all 
Xi,.... X rn and all matrices A which satisfy the assumptions, the convergence in (4.45) follows. The result 
may thus be regarded as a multivariate version of the direct half of Karamata’s Theorem. 

Proof of Proposition 4.4. Define independent random variables XI, 1 < j < m, such that XI has P x i- 

density x l[o,oo )( x ) x ^ (E(Xj 3 J ) -1 ,1 < j < m. This is possible because all /3j G [0,1) and thus E(Xj 3 ) < 
oo by our assumptions. For those 1 < j < m with Kj > 0 the random variable XI is regularly varying with 
index — (1 — ft j), because 


P(X'>x) f x °°y*P x *(dy) 

™ xPiP{Xj > x) E{XPi)xhP{Xj > X) 


= (E(X*)(1-I3 j ))- 


for all 1 < j < m by Karamata’s Theorem (cf. Cline (1983), Lemma 1.1). Thus, for 1 < j < m with Kj > 0 
the random variable Xj := (X') 1_ ^> i s regularly varying with index —1 and 


lim P (Xi> xl fe ) = (E{Xb)( 1 - fa))- 1 , 
x^oo xhP(Xj > x) y y A 


(4.46) 


For 1 < j < m with Kj = 0 we have Xj > 1 a.s. because we assumed Xj > 1 a.s. Furthermore, we have for 
all S G (0,1] and all j with Kj = 0 that 

dx) _ E{X)- (1 -^ 


E{X)~ S ) = 


< oo. 


E(X”) E{XV) 

Thus, the random variables Xj , 1 < j < m, satisfy the assumptions of Proposition 4.1. Set now 

A =(<*«):= ((l-ftrVi)GR" Xm 

Then k := ((1 — 0j)Kj)i<j<m is the unique solution to the linear program 


find x > 0 such that Ax >!. E 
and Ak = 1. Set 


2—1 


Af(A) = { (an 


, . . . , Xm) : 

1=1 


Then, 




Xj' 3 > 1,1 < i < n 


(Xj )i <j<m G y d/(A) 

>1, 1 < i < n, 


ni 


1=1 


y" 
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O- 


and so 


m 


n 


1 y(l~fij) K j 


(1 /3 j) dij 

>1, 1 < i < n, 
(Xj)l<j<m £ V ®n M(A), 


'([Xj)l<j<m€y®KM( A)) = J p( X j)l<i<™( dx) 


y0/*M( A) 


ji 


J 

y®*M(A) J J 

f TT (2/ Kj p(Xj/y"J )i<i<m ( dx ). 

J E(X^) 


M( A) 


Thus, 


M(A) 

lim -—- 7 --- 

y->°o P(Xj > y K i) 

i-Kj> o 


= lim 

y—too 


P((Xj)i<j<m £y® A M(A)^ 

n P{X j>y *s) 


n^(xf) n p{x j>y *i) 

j =1 j-Kj> 0 


m ’ 

II //"' v - 11 P{Xj > y Kj ) 

3 =1 j:Kj>0 


where the first factor converges to 


/ 

M(A) 


II x d 2 \{d(xj) {j .. K . >0} ) <g> p(^ ) 0' : ^- 0 >(d(^) {j:K . = o } ) 

j:Kj>0 


(4.47) 


(4.48) 


by the assumption. Substitute (z/i,..., y m ) := (a^ 1 ^ ,..., a;^ /3m " ) ) and note that (xi,..., x m ) € M(A) 

is equivalent to (j/i,..., y m ) € M(A), so the expression in (4.48) equals 


/ II 2 A(d( % ) b - :Kj> 0 } )®p ( ^ ) f^ J =o, (d(%){j:Kj=o}) 

M(A) ■^ >0 

= / f n ( n 

M(A) V^ >0 / V :K *=° / 

A ( d (yj){j:^>0}) ® P ( ^ ) «-J'=0}(d(2/ : ,) {j:K;) . =0} ). 

The second factor in (4.47) converges to 

ru =t , E(*f) 

— fij) 

by (4.46). Taken together, this yields the statement of the proposition. □ 
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